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Abstract 

An experimental and theoretical study on the structure and energetics of stellated gold clusters 
at several sizes is presented. Systematic molecular dynamics simulations on Kepler-Poisont classi- 
fied clusters are performed based on cubo-octahedral and icosahedral cores, following present and 
previous studies that suggest stellated clusters grew up from these seeds. For cluster sizes up to 
10802 atoms full atomistic molecular dynamics simulations at room temperature have been carried 
out. Results show that stellated clusters in the bigger size regime maintain their star-like shape at 
room temperature, and by means of linear fitting of energy data to proper function models, it is 
predicted that the two type of stellated structures might coexist at room temperature for relatively 
large sizes, being these findings in good agreement with experimental results. 

PACS numbers: * 



I. INTRODUCTION 

The present study on stellated gold nanoparticles (NPs) is driven mainly by fundamental 
and applied leit motivs. From the fundamental point of view, we pursue the hypothesis that 
stellation on clusters (Kepler-Poisont solids) is a process that follows after the nucleation of 
well defined structures, such as octahedra and icosahedra (platonic shapes) and cuboctahe- 
dra and truncated octahedra (archimedean solids) 1 , supported the idea by the seed growth 
approach in many different syntheses of anysotropic noble metal NPs 2-7 ; where each one of 
these shape classifications could be located in size regions — althought with not well defined 
boundaries — of the cluster growth. 

From the applied point of view, in the nanoparticle research fields there are now well stab- 
lished facts that drive applications researchs in these fields, i.e., the properties of the systems 
are influenced by the NPs size, but also is known now that the shape strongly influences 
the physical and chemical properties of the NPs 8-10 . Thus, in the problem of understanding 
how shape and size of metallic NPs relate each other, and how these factors influence the 
properties of the systems under study, becomes more complex but certainly, more fascinating 
too. Because of the facile synthesis, their enhanced absorption and scattering cross sections 
that Au and Ag NPs show under electromagnetic radiation, and because of their poten- 
tial nontoxicity, noncitotoxicity, and their good biocompatibility, they are widely studied 
for biological imaging and biomedicine applications (cancer cell imaging and destruction, 
photermal therapy, among others) n ~ 13 . In particular, Au NPs generate heat under proper 
electromagnetic radiation, being this heating effect strong when the energy of the incident 
radiation is close to the plasmon frequency of the metallic NPs 14 . 

Among factors that had contributed to the comprehension, understanding and develop- 
ment of the energetics, thermodynamics, and kinetics of atomic and molecular clusters 15 , 
we can mention the interplay between spectroscopy 2 , and theoretical techniques based on 
solid-state physics (SSP) 16 ~ 19 or other potentials supported by the density-functional the- 
ory (DFT) 20 ' 21 . The development of geometric models based on platonic solids 22 , or on the 
most probably structures in experiments generated by quenching technique in a molecular 
dynamics simulation 23 or by genetic algorithms 24 provide us with suitable geometrical and 
inter-atomic potential models are available to study small and large number of particles 
assemblies through non-quantum methodologies like Monte Carlo or molecular dynamics 



simulations (MDS) at high confident levels. 

Under this theoretical framework, extensive studies on metallic clusters have been carried 
out 15 . For instance, by using the excess-energy per number of surface atoms on clusters as 
a diagnostic tool, Balletto et al. 25 had found and explained why gold clusters prefer fee 
structures at sizes larger than 600 atoms, and strong competitions between D^ and fee near 
N = 400 atoms. Two potential models based on DFT and SSP theories were used in those 
calculations, and they also found agreement with other results obtained at different levels 
of calculation. 

Several potential models for transition and noble metals, based on the tight-binding 
second-moment approximation (TBSMA) 16 ' 17 , have been proposed. For instance, the pa- 
rameterization of Rosato et al. ls (RGL model) and a second parameterization of the RGL 
due to Cleri and Rosato 19 (CR), among others 25,26 . All these kind of models include the 
repulsive pairwise Born-Mayer potential to stabilize the system. One of the main charac- 
teristics of those models is the inclusion of many body effects ensuring correct estimation 
of several bulk properties, like cohesive energies, vacancies, surface effects, etc. The correct 
estimation of those bulk properties, mainly cohesive energy, surface effects, and surface re- 
costructions, has given confidence to extend the model to atoms in a cluster environment 15 . 
Thus, with these reliable potential models at hand, general trends can be obtained. For 
example, significative knowledge on the thermal behavior of gold clusters is well known for 
specific geometries based on crystalline and non-crystalline motifs, like fee, icosahedral and 
decahedral clusters 15 . 

In this work, by means of classical MD simulations, and the CR interatomic potential 
model, we study the structure and dynamics of stellated gold nanoclusters with icosahedral 
and cuboctahedral seeds, biased mainly by the present experimental results on gold nanopar- 
ticles and previous reports on stellated (pinata like) clusters 2 . Stellated Au nanoparticles 
have been synthesized at room conditions of temperature and pressure. Figure 1 shows the 
great rate of production of these stellated clusters, and the high-resolution transmission elec- 
tron microscopy (HRTEM) images at different orientations clearly show highly symmetric 
pyramidal structures with slightly smoothed apices. From experimental facts, it has been 
observed high stability in the clusters shapes, i.e., they do not change in days and even 
months after depositing them in an inert environment. These facts induced the hypothesis 
that stellated gold particles could have cuboctahedral and icosahedral cores, but both with 



a surface tetrahedral pyramidal growing along the high symmetry axes. In order to give 
some insight into the structural and thermal behavior of these stellated structures, mainly 
at room temperature, we present a systematic study of the structural and thermal behavior 
of two groups of gold clusters whose construction was based on those structures. 

II. EXPERIMENTAL DETAILS 

A. Materials and Methods 

Hydrogen tetrachloroaurate trihydrate (HAuC^ • 3H2O), Silver nitrate, trisodium citrate 
dehydrate, ascorbic acid were purchased from Sigma-Aldrich and used as received. Multi- 
branched polyhedral gold nanoparticles were prepared by colloidal reduction in aqueous 
solution at ambient conditions according to previous literature methods 2 . A volume of 0.5 
mL of 0.1 M HAuCl 4 aqueous solution (50 /xmol) was rapidly added into an aqueous solution 
of ascorbic acid (2 mmol in 40 mL H 2 0) under vigorous stirring. The solution immediately 
changed its color to opaque orange-red within 1^2 seconds. All the samples were purified 
by washing with ethanol and centrifuging, and finally dispersed in ethanol. 

B. TEM Characterization 

The Au nanostar samples were characterized by ultra-high resolution scanning electron 
microscope (SEM) FEG Hitachi S-5500 (0.4 nm at 30 kV) with BF/DF Duo-STEM detector 
and high resolution transmission electron microscope (HRTEM) Jeol JEM-2010F FE-TEM 
with an accelerating voltage of 200 kV, equipped with a Schottky-type field emission gun 
and a point-to-point resolution of 0.19 nm. Samples for TEM studies were prepared by 
placing directly the gold nanoparticles onto a holey carbon TEM grid. 

III. THEORETICAL BACKGROUND AND METHODOLOGY 

The functional form of the TBSMA for pure atomic systems, also known as Gupta 
potential 16 , is given by 

V(n) = ^Aexp(-pAry) - I ^£ 2 exp(-2g Ar 4i ) J . (1) 



The first term is the repulsive Born-Mayer pairwise potential, and the second one is the tight- 
binding based attractive non-pairwise potential which takes into account many-body effects. 
Ar^ = {rij — r )/r is the fractional change in distance between atom i and atom j, and r is 
the distance between first nearest-neighbors in bulk. The five parameters, A, £,p, q, and tq 
must be determined for each particular transition or noble metal atom by a parametrization 
procedure using bulk data 18 ' 19,27 . If only first nearest-neighbors are taking into account in 
Eq. (1), and if the equilibrium condition is fulfilled, then the number of parameters is 
decreasing from five to three 17 : 

A = -*-£?*, (2a) 

p — q Z 

£=— ^, (2b) 

p-q Vz 

where _E co h is the experimental cohesive energy of bulk used in the parameterization, and Z 
is the number of first nearest-neighbors in the crystal. 

Some sets of parameters for gold (fee crystal) are given in Table I. In these models values 
for A and £ are very close to those calculated from Eqs. (2) using the experimental value 
E co h = —3.78 eV/atom and Z = 12, both from gold fcc-crystal 28 . At a first sight, differences 
and similarities between these potential models can be obtained from the fractional change 
in total energy accompanying a fractional change in distance between first nearest-neighbors 
e = ( r ij ~ r o)/ r cb which at second order approximation in e is roughly given by 

AE 



\E(0)\ aie ' (3) 

where AE = E(e) — E(0), E is the total energy, and o\ = pq/2, see Ref. 25. The o\ 
values for the CR and RGL models differ 1.5% each other, and the values among BFFMM 
and CP differ 0.8%; corresponding the lower and bigger values to the CR and CP models, 
respectively. Thus, according to 0\ values and Eq. (3), the BFFMM potential is as sticky 
as the CP model, and the CR model is slightly less sticky than the RGL, being CP more 
sticky than the CR model. However, according to the formula 

A E 1 

-(q(e-^-l)-p(e-^-l)), (4) 



\E(0)\ p-q 

where no approximations in e have been used, the similarity between the BFFMM and CP 
models only occurs for displacements less than 1.5% of r^, as can be seen from the plots in 
Fig. 3. In fact, crossovers between the CP and the BFFMM, RGL and CR models occur at 



TABLE I: Comparision of published potential parameters for gold derived by fitting bulk experi- 
mental data to values calculated using the Gupta functional form given by Eq. (1). We note that 
values for the lattice constant and cohesive energy, (y/2 ro, -Ecoh)) i n each model are experimental 
values taken from Kittel 28 , but on different potential versions: the parameterization of Rosato et 
al. 18 (RGL model, Kittel 1972); the second parameterization of the RGL due to Cleri and Rosato 19 
(CR, Kittel 1966); whereas Balleto et al. 25 (BFFMM) and Chamati and Papanicolaou 26 (CP) use 
parameters from Kittel, 1996. For the CP model we have used the calculated values 26 shown in the 



row. 



Model a 


P 


q 


£(eV) 


A (eV) 


V2 r (A) 


-E coh (eV/atom) 


RGL 18 


10.15 


4.13 


1.8398 6 


0.2161 6 


4.08 


3.78 


CR 19 


10.229 


4.036 


1.790 


0.2061 


4.079 


3.779 


BFFMM 25 


10.53 


4.30 


1.855 


0.2197 


4.07 


3.78 


CP 26 


14.6027 


3.12572 


1.32795 


0.08170 


4.068 


3.76 



"Model's name is proposed from the initials of the author and coauthors of the cited article. 
^Values obtained using Eqs. (2). 

e ~ 0.0085, e ~ 0.093, and e ~ 0.115 respectively. Thus, according to Eq. (4) the CP model 
is less sticky than the other potentials, but for e > 0.115. 

A. Simulation Details 



In order to perform MD simulations at room temperature of the SC and SI gold clus- 
ters described above, the Au potential model given by Eq. (1) was used, applying the four 
parameter sets given in Table I under the GROMOS96 scheme 29 . This package applies the 
leap-frog algorithm to integrate the Newton's classical equations of motion under vacuum 
or periodic boundary conditions 30 . High numerical stability is ensured for a time-step of 2 
fs, if shake method is applied; otherwise 0.5 fs. Those time-steps values are commonly used 
for flexible biological molecules 29 ; for metallic atoms a practical value for the time-step is 1 
fs 31 . Also under this scheme, the system could be weakly or strongly coupled to a temper- 
ature bath using the Berendsen's thermostat 32 . The numerical ensemble produced by MDS 
applications using Berendsen's thermostat is the weak-ensemble 33 , which has intermediate 



property between the canonical (JV, V, T;a = 0) and microcanonical (N,V,E;a = 1) en- 
sembles. Here, a is a constant in equilibrium but dependent on the temperature coupling 
parameter 32 ' 33 . In the practice, a is determined from the ratio of standard deviations of 
kinetic and potential energies, 



a 



^((5E Km ) 2 )/((5E Pot ) 2 ), (5) 



where 5E z (t) = E z {t) — (E z ) stands for the instantaneous fluctuation of property z (kinetic 
or potential energy), and ( ) denotes average time. 

B. Geometrical Description 

Since any stellated gold clusters considered in this work are made up of one icosahedral 
or one cuboctahedral cluster in the core, and regular tetrahedral clusters attached to their 
(111) faces, we will refer them as stellated icosahedral (SI) and stellated cuboctahedral (SC) 
clusters, respectively. Specifically, each SI gold cluster is made up of one Mackay icosahedral 
cluster with {n — k + 1, k ^ 1) shells and N(n) = yri 3 — 5n 2 + yn — 1 atoms, and twenty 

tetrahedral clusters, each with N<?d(k) = — g atoms and attached to each one of the 

twenty (111) facets of the icosahedral cluster. In a similar way, each SC gold cluster is made 
up of one cuboctahedral cluster with N(k+1) atoms, and eight tetrahedral clusters, each with 
iVrd(fc) atoms and attached to each one of the eight (111) facets of the cuboctahedral cluster. 
Hence, the number of atoms of any SI^ cluster is given by N$i(k) = N(k + 1) + 20N^^(k), 
while N$c{k) = N(k + 1) + 8iVTd(&) is the number of atoms of any SC^ cluster. N^i(k) and 
Nsc{k) are only different in the number of tetrahedrons, and they can also be expressed as 

14 19 

Nsc(k) = — k 3 + 9k 2 + Y k + 1 > ( 6a ) 

on ^1 

N sl (k) = -k 3 + 15k 2 + -k + l, (6b) 

with k ^ 1. 

The gold clusters considered in this work are those SC clusters with k ^ 13, and those 
SI with k ^ 11, shown in Fig. 4. In all the initial geometrical models the distance be- 
tween first nearest-neighbor atoms is fixed to the first neares-neighbor distance of bulk gold, 

r = 2.88 A. 



IV. RESULTS AND DISCUSSION 

Figure 4 displays the unrelaxed stellated clusters used in this study. To prepare these 
SI and SC gold clusters for MD simulations at T = 298 K, we have first slightly relaxed 
each one of them with the steepest-descents method using the four potential models in 
Table I. Although the main purpose of the minimization procedure is to reduce the internal 
strain of the cluster, and the consequently reduction of computer effort at the equilibration 
period, we have analyzed these results to gain insight investigating how the residual stress is 
reflected into the structures, how the surface and volume contributions to the total energy 
behaves as function of cluster size, and what type of the stellated structures is the most 
favorable. Similarities and differences between the results obtained using older and newer 
gold potential models (Table I), are also enhanced and analyzed at this stage in order to 
choose the appropriate potential to perform the simulations at room temperature. Thus, a 
conciencious analysis of the structural and energetic results obtained under this procedure 
is presented in Sec. IV A. In Sec. IV B, we have posed similar questions for the equilibrated 
clusters at T = 298 K, and the results have been analyzed in the same way, but only for the 
CR potential model, which was elected after the analysis in Sec. IV A. 

A. Relaxed Clusters 

In cluster physics is common to describe atoms as hard spheres, and to use them as a 
linear lengt unit. Following this idea, in the micro-crystal size regime (where a high rate of 
the total N atoms belong to the bulk) the particles are iV 1 / 3 atoms long, with area facets 
proportional to N 2 ^ 3 atoms. Only some atoms belong to the surface (~ iV 2 / 3 ), located 
predominately on the facets and only a negligible number on the edges and the vertices. 
On the contrary, when the cluster size is decreased, almost all the atoms tend to be at the 
surface. In this case, however, the amounts of atoms in the edges, the vertices, and the 
facets tend to compete among them, and the extreme case is for the very smaller clusters 
that contains only vertex atoms. Therefore, it would be helpful to obtain a function that 
give insights in which range of N, the macro/micro (iV sur f <C N) or the nano (iV 6ur f ~ N) 
character of any studied cluster is manifested. 

For this purpose, it is proposed a function for the number of surface atoms iV sur f, in terms 



TABLE II: The parameters for the fitting function, Eq. (7), together with its asymptotic standard 
errors for the SC and SI gold clusters relaxed with the steepest-descents and the CR potential 
model. 

a (3 7 

SC 6.4490 ± 0.0005 -13.97 ±0.01 9.33 ± 0.08 

SI 8.4663 ± 0.0006 -39.74 ±0.02 73.15 ±0.11 

of iV 1//3 , iV 2 / 3 , and a constant, as follows 

N smi (N) = aiV 2/3 + /3iV 1/3 + 7 , (7) 

where the coefficients a, (3, and 7 must be free parameter, since the function must describe 
the full number of surface atoms for small, intermediate, and large size clusters. Thus, when 
fitting the function to data a high correlation between them is expected and some of them 
could get negative values. 

The behavior of the computed ratio, N SUT f/N 2 ^ 3 , as a function of the cluster size, is almost 
the same for both the relaxed and unrelaxed SI and SC clusters, as is shown in Fig. 5. The 
number of surface atoms of each relaxed cluster, N sm {, has been calculated through the Lee- 
Richard algorithm 35 using the Van der Waals radius of gold as a probe radius r probe = 1.66 A. 
Hereafter we denote this function as J\f(N) = iV sur f/iV 2//3 in order to facilitate discussions. 
In the Figure, the solid and dotted lines are obtained by linear fitting of SC and SI N sur f 
data to functions of the form given by Eq. (7), respectively. Here the function J\fsc{N) is 
for iV ^ 21, while Afsi(N) is for N ^ 135. The coefficients of the fitting functions for the SC 
and SI clusters are given in Table II. Note that some parameters are necessarily negative. 

As can be seen from Fig. 5, calculated values for SC clusters tend faster to its asymptotic 
value, ctsc = ±6.449, than those values for SI clusters, that tend to «si = ±8.466. The 
asymptotic value «sc is greater than the asymptotic value for a simple cubic cluster, a cu bi C = 
6.0, since the SC clusters are based on a truncated fee with tetrahedral units attached to 
the (111) faces. The value «si = ±8.466 is explained by a similar reasoning. For cluster 
sizes of O(10 6 ), values of the fitting functions A/ sc and A/si are around the 2.2% and 4.6% of 
their asymptotic values, respectively. Thus, one could expect stronger surface effects for SI 
clusters than for the SC clusters, or one could expect that the SC clusters already behave 
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as macro clusters (iV sur f oc iV 2 / 3 ) for relatively small sizes, N ~ 10 6 , so that many properties 
of SC could be dominated by volume atoms. For SI clusters in the same range one could 
expect properties dominated by surface atoms. As a note of caution, the stellated clusters 
have very specific sizes, see Eqs. (6), because they are constructed using complete tetrahedral 
units mounted on the eight triangular faces of the cuboctahedra or over the twenty faces of 
the icosahedral clusters. Hence, detailed interpolation or extrapolations results from Eq. (7) 
must be taken with care. 

In the small size range, the stellated clusters with sizes iVsc(4) = 469 and iVgi(4) = 709 
contain tetrahedral units of 20 atoms, and interesting is to note that tetrahedral gold cluster 
with 20 atoms have been proven to be the more stable gold isomer 36 . To this respect, we 
also note that after the slightly relaxation processes, where each of the potential models 
displayed in Table I were used, we have observed that the shape of each cluster in Fig. 4 
is maintained; the relaxed clusters are contracted in the average, as it is expected for these 
kind of potential models 16 , being more evident this result in the smaller ones due to the high 
value of the surface /volume ratio, but the peaks of all the clusters are maintained, which 
can be attributed to the same mechanism that maintains vertex atoms of the global minima 
tetrahedral cluster of Fernandez et a/. 36 . 

In order to explore the energetic trends of the relaxed SC and SI clusters under the 
potential models given in Table I, we have first fitted the potential energy data of each 
series to a function of the form (see Ref. 15 and references therein) 

E(N) = AN + BN 2/3 + CN 1/3 + D, (8) 

where the coefficient A is proportional to the energy per atom of the bulk-like atoms, that 
follows from the asymptotic behavior of E(N)/N for large clusters. The number of bulk 
atoms is taken proportional to the cluster size N. Thus, at least for large clusters, the first 
term in the Equation could be interpreted as the contribution to total energy due to bulk 
atoms, and the last three term together as the contribution to total energy due to the atoms 
at the surface. The explicit dependence on N n (1/3, 2/3, and 0) comes from the splitting 
of the full number of surface atoms into facets atoms (proportional to iV 2 / 3 ), edges atoms 
(proportional to iV 1 / 3 ), and vertexes atoms, a constant. 

Therefore, by construction, the first term in Eq. (8) is interpreted as the total volume 
contribution to the energy, and the second, third, and fourth terms together are interpreted 

11 



TABLE III: The coefficients for the fitting function, Eq. (8), for the steepest-descents relaxed SC 
(first row) and SI (second row) clusters with the potential models given in Table I. Units are 
eV/atom. 



Model 


A 


u a 


B 


C 


D 


CR 


-3.7795 ± 0.0003 


0.0000 


1.570 ±0.007 


-0.71 ±0.06 


1.68 ±0.11 




-3.7687 ±0.0002 


0.0108 


1.678 ±0.003 


- 


-3.53 ±0.05 


RGL 


-3.8256 ±0.0002 


0.0000 


1.454 ±0.007 


-0.73 ±0.05 


1.76 ±0.10 




-3.8160 ±0.0002 


0.0096 


1.554 ±0.003 


- 


-3.21 ±0.04 


BFFMM 


-3.8284 ±0.0002 


0.0000 


1.436 ±0.005 


-0.77 ±0.04 


1.84 ±0.08 




-3.8186 ±0.0002 


0.0099 


1.533 ±0.004 


- 


-3.16 ±0.05 


CP 


-3.7600 ± 0.0001 


0.0000 


2.967 ±0.004 


-0.54 ±0.03 


0.34 ±0.06 




-3.7336 ±0.0013 


0.0264 


3.125 ±0.021 


- 


-7.87 ±0.29 



V = A - E Tei , where E Tei = A sc . 

as the total surface energy contribution (see Ref. 15 for details). As in the case of N sm f(N), 
see Eq. (7), the coefficients in Eq. (8) must be free parameters, and high correlation among 
them are expected (mainly among coefficients B, C, and D) when fitting data of small, 
intermediate or large size clusters to this equation. 

The parameters for the fit to each type of stellated clusters and potential models are 
given in Table III. To avoid large asymptotic standard errors in the linear fitting data of SI 
clusters, the coefficient C is excluded from its fitting process. Note that, for the potential 
models RGL and BFFMM, the asymptotic values for the binding energy, Asc, are lower 
than the experimental cohesive energy of bulk used in the parameterization of each of the 
potential model (see Table I), however they are very close to the experimental cohesive 
energy, E co ^ = —3.81 eV/atom reported in Ref. 37. For the CP model Asc is exactly the 
calculated value of the cohesive energy of bulk Au reported by Chamati and Papanicolaou 26 , 
and for the CR model the value of Asc is very close to the experimental cohesive energy of 
bulk used in its parameterization 19 . 

In order to explore which type of relaxed structure is more favorable, we have used as 
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diagnostic tool the excess energy per iV 2//3 , defined by 15 , 

. / E — NE rp r , . 

A ( N ) = ^2/3 - ; ( 9 ) 

where E stands for the total potential energy of one cluster of size N, and it is common to 
take for the reference energy E Ie f the experimental cohesive energy of the bulk used in the 
parameterization of the potential model, see Table I. Here, for each potential model, A(iV) is 
calculated using E rei = A S c, where A sc is given in Table III. The results of the calculations 
of A (AT) using each of the potential models in Table I and the steepest-descents method 
are given in Fig. 6. It can be observed from this figure that st this stage of the simulations, 
almost all SC clusters, except iVsc(l) — 21, are more favorable than the SI clusters. This 
behavior occurs for all potential models, except the CP model which also excludes iVsc(2). 
Similar results have been found by Baletto et al. 25 , using the RGL and embedded-atom 
models, but for non-stellated icosahedral and cuboctahedral clusters. However, as can be 
seen from the calculated data in the figure, the effects of mounting tetrahedral motifs on the 
triangular faces of cuboctahedral or icosahedral clusters changes significantly the behavior 
of A(iV). For instance, A data of icosahedral clusters in Ref. 25 is concave up, while the A 
data of SI clusters in Fig. 6 is concave down. For Lennard- Jones icosahedral and fee cluster 
data 38 ' 39 , it has been found that A is concave down for both types of clusters. We have 
explored the concavity effect by expressing A(iV) in the form 

A(JV) = aN l/3 + B + C/N l/3 + D/N 2/3 , (10) 

which follows using Eq. (8) into Eq. (9). The parameter a = A — E Te f is interpreted as the 
strain or elastic energy per atom of the cluster 15 ' 39,40 . It is well known that for clusters with 
fee symmetry, <jf cc is negligible since the bulk energy per atom A{ cc , is the optimum. On 
the contrary, it is also well known that the non-crystalographic symmetry of the clusters 
frustrates the elimination of strains by the minimization process, and consequently they can 
accumulate that energy as elastic-energy 15 ' 39 ' 40 , which increases with cluster size. In principle, 
A must be higher or equal to the cohesive energy of the bulk crystal, -E C oh, which can be 
appreciated as the optimum energy per atom of volume atoms; however in practice, when 
fitting MDS results to Eq. (8), A could be slightly lower than -E co h- This could be attributed 
to the linear fitting of MDS data, and or to the cut-off radius used in the interactions: in all 
potential models we have included more than one coordination shell in the interactions. In 
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any case, since one is interested in the relative stability of the two types of clusters, to avoid 
negative values for a in Eq. (10) we have used for the reference energy E rc f the lower value 
for A in each potential model, instead of the experimental cohesive energy of bulk used in 
the parameterization of the potential model, see Table. I. 

In Fig. 6, the solid and dotted lines in each plot are the curves Asc(^) and Asi(N) 
obtained from Eq. (10) with the coefficients given in Table III, for the potential models 
CR (a), CP (b), RGL (c), and BFFMM (d). Note that, in general, the asymptotic value 
for the binding energy Age, is lower than the experimental cohesive energy of bulk used in 
the parameterization of the potential model, see Table I. As pointed out in the preceding 
paragraph, to avoid negative values of crsc, we have used E Te f = A$c in the Eq. (10) for the 
SC and SI clusters. This procedure has not effect in the present analysis, since only we have 
changed the reference point of energy. For all the potential models the respective functions, 
Agc(iV) and Asi(iV), predict only one crossover size. Only the crossover sizes for the RGL, 
BFFMM, and CR models agree with the rule: the larger is pq, the smaller are the crossover 
sizes 15 ' 25 . The CP model has the higher value of pq, and also the higher value of the crossover 
size. The exact values of the crossover sizes have no real significance, since the calculations 
have been limited to clusters sizes given in Eqs. 6. We only could infer from these results 
that both types of clusters could coexists around these sizes. For any potential model, no 
others crossover sizes occur for the larger clusters of both types, since their asymptotes, right 
lines in A(iV 1//3 ), are not intersected. Thus, each one of the potential models predict only 
one crossover size, which occurs for relatively small cluster sizes. Therefore, at this stage of 
the MD simulations we could partially conclude that (a) all the slightly relaxed structures 
maintain the star-like shape of the initial structures; (b) all SC clusters, except the smaller 
cluster, are more favorable than the SI clusters; (c) the most strained SI structures are given 
by the CP model, and the less strained are given by the RGL model (elastic energy values 
given by the CR and BFFMM models are very close to the value given by the RGL model); 
and (d) only the asymptotic value for the binding energy Asc, given by the CR model, is 
very close to the value of the experimental cohesive energy used in its parameterization. 
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B. Equilibrated Clusters at T = 298 K 

MD simulations at room temperature of the relaxed clusters described above have been 
conducted as follows: for each cluster a sample of mechanical states of length 0.25 ns has 
been produced in a weak-ensemble approximating the microcanonical-ensemble (a ~ 1), 
and after 2 ns of equilibration length. This procedure has been applied to each cluster in 
the range 21 ^ iV ^ 3095, including N = 5631 and N = 6169, and except those with 
sizes iVgc = 9437 and iVgi = 10803, where an equilibration length of 1 ns was used. In the 
equilibration period, the Berendsen's thermostat with time-coupling parameter tt = 0.4 
ps and target temperature T = 298 K was used 32 ' 33 . In all the simulations a time-step of 
1 fs is used, and periodic boundary conditions has been applied, but with a much larger 
box length than the cluster radius to ensure non-interacting images. A cutoff radius equal 
to half box length was used for clusters sizes N ^ 87, otherwise a cutoff radius of 2 nm 
was set. Statistical properties of the produced sampled distribution were obtained using the 
coarse-graining stratified systematic sampling dividing the sample in equal bins of length 
0.05 ps, to avoid serial correlations. The averages of the properties considered in this work 
are reported at the 95% of confident level. All simulation results have been obtained using 
only the CR potential model, since the asymptotic behavior of E(N)/N for the relaxed SC 
clusters reproduces very well the experimental cohesive energy used in the parameterization 
of the model. 

As it can be observed from the average equilibrium configurations of the two series of 
clusters shown in Fig. 7, one of the more evident temperature effects on clusters for sizes 
347 < N < 9437 is the smoothing of all the peaks of the tetrahedral pyramidal motifs of the 
relaxed clusters. This is easy to understand due to the natural tendency of lower coordinated 
surface atoms to increase its coordination number to the optimum value, which in this case 
is driven by thermal energy. For the first three clusters of each series we have found that 
all the peaks are completely destroyed in agreement with recent results, which shown that 
small gold clusters are low-symmetry or amorphous-like clusters 40 . By the other hand, for 
the two largest clusters of each series, Figs. 11, some peaks remain due mainly to the length 
of 1 ns used in its equilibration period. Their contribution to the energy could be neglected 
since, as we will show below, volume contribution dominates over surface contributions at 
these sizes. Eventually, it is expected these peaks are finally smoothed if equilibration length 
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TABLE IV: Same as Table II, but for the equilibrated gold clusters with the CR potential model. 

a f3 7 

SC 6.41 ±0.02 -13.07 ±0.33 

SI 8.44 ±0.03 -37.32 ±0.50 

is extended to 2 ns. Thus, the surfaces of our equilibrated clusters at room temperature for 
sizes in the range 347 < N ^ 10803 resembles the TEM images of the gold particles (a), (b) 
and (d) in Fig. 2. 

By following the equilibration time history of clusters with sizes greater than 347 atoms, 
we have observed in detail how the peaks of the pyramids are smoothed. For any cluster 
at the initial steps of the equilibration process, each pyramid tip atom is diffused along one 
edge and it attaches around the middle edge for long time, generating metastable states. 
After a short period, these atoms remove (almost always) another atoms from those mid 
edge positions, that are rolled towards the base of the pyramid, and the total energy is 
stabilized for longer periods, which might be considered as an stable state. For the SC 
clusters, we have found that those atoms form oscillating pairs of atoms, which are deposited 
on the surface around vertices connecting square faces. Nowdays, there are high interest on 
the collective electronic excitations or surface plasmons of star-like gold nanoparticles, and 
these theoretical findings could give some insights on the features of the collective electronic 
excitations resulting from the interaction of star-like nanoparticles with the electric field of 
light 41 " 44 . 

From the above results on equilibrated surface atoms, minimal changes in the behavior of 
the number of surface atoms as a function of cluster size are expected, and this is indeed the 
case if the curves J\f(N) for the SC and SI clusters on Fig. 5 are compared with those of the 
equilibrated clusters shown in Fig. 8. The behavior of the corresponding curves, for N > 347, 
is very similar. As in the case of the relaxed clusters described in the preceding section, curves 
in Fig. 8 have been obtained by linear fitting N sut { to Eq. (7), where iV 8ur f is the number of 
surface atoms of the average equilibrium configurations shown in Fig. 7. The coefficients of 
the fitting functions for the equilibrated SC and SI clusters are given in Table IV, where the 
coefficient 7 has been excluded from the fitting process in order to avoid large asymptotic 
standard errors in the coefficients. Since the main interest here is on the large stellated 
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TABLE V: Same as Table III, but for the equilibrated gold clusters with the CR potential model. 



Cluster type 


A 


a a 


B 


C 


D 


SC 
SI 


-3.6982 ± 0.0006 
-3.7132 ±0.0040 


0.0151 
0.0000 


1.496 ±0.012 
2.545 ±0.107 


-8.22 ±0.63 


-1.64 ±0.54 



V = A — E rc {, where E rc { = A 



si- 



structures which maintain their peaks, the first three equilibrium configurations of each 
series have been excluded from the fitting, and also from the discussions, since they are 
amorphous or low-symmetry structures described in detail by Soler et al. 40 . ^From Fig. 8, 
the number of surface atoms of the SI clusters is lower than for relaxed clusters; on the 
contrary, relaxed and equilibrated SC clusters have almost the same number of surface 
atoms. For large clusters, the asymptotic behavior of dotted curves, or solid curves, is very 
similar. Thus, from the energetic point of view, details of the surface atoms of very large 
particles could not be important since volume atoms dominate energy contributions, but the 
smoothed peaks of the clusters can play a significant roll in the description of the surface 
plasmons of these nanoparticles. 

By following the same procedure as for the relaxed clusters, we have also explored which 
type of equilibrated gold cluster is more favorable. In this case E in Eq. (9) is the total 
energy of a cluster of size N, N suri is given by Eq. (7) and Table IV, and for E rci we 
have used A$i = —3.7132 eV/atm given in Table V. The calculated values for Asc(N) and 
Asi(iV) are given in Fig. (9) together with its fitted curves, where the coefficients are given 
in Table V. In this case only E(N) data for iV ^ 347 has been fitted to Eq. (8), and to avoid 
large asymptotic errors in the coefficients, C and D have been omitted in the fitting process 
of data E SC (N) and E S i(N), respectively. ^From interpolation and extrapolation results we 
have found that the SC clusters are more favorable in the range 671 < N < 224 586, while SI 
clusters are more favorable out of this range, excluding the first five smaller clusters. Thus, 
we predict that stellated gold clusters coexist around two crossover sizes; one at N ~ 671 
and another one at N ~ 224 586, a result that is well agreement with experimental facts 
(see Fig. 2). 
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V. CONCLUSIONS 

We have synthesized and characterized by TEM stellated gold nanoparticles, as well as 
performed a systematic study for their structure and energetic stability through geometrical 
models and by means of molecular dynamics simulations at room temperature. From the 
analysis of the MDS results, we can conclude that: (a) the morphology and the surface struc- 
ture of these stellated clusters for sizes N > 347 are in good agreement with the synthesized 
NPs, as well as with the results reported recently by Burt et al. 2 . The structure of peaks is 
maintained and the coexistence of the two type of particles has been demonstrated. Smaller 
clusters are low-symmetry or amorphous structures already studied in Ref. 40; (b) the stel- 
lated gold clusters prefer SC tructures for cluster sizes in the range 347 < N < 226281 
(« 30 nm), where surface effects are stronger for SI clusters than for SC clusters; in fact the 
behavior of Asc(^ 1//3 ) is almost lineal. We predict that around a cluster size of 30 nm, both 
types of clusters could coexist, and for larger clusters the SI structure is more favorable than 
the SC structure. Finally, it is important to emphasize that we are comparing clusters of 
different size which have cores of the same size, cuboctahedral or icosahedral clusters with 
same magic size, thus this difference in the number of atoms arises from the tetrahedral 
motifs attached to their triangular faces. 

The present study on stellated gold nanoparticles is driven mainly by fundamental and 
applied leit motivs. From the fundamental point of view, we advanced the hypothesis that 
stellation on clusters (Kepler-Poisont solids) is a process that follows after the nucleation 
of well defined structures, such as octahedra and icosahedra (platonic shapes) and cuboc- 
tahedra and truncated octahedra (archimedean solids), where each one of these shape clas- 
sifications could be (althought with not well defined boundaries) located in size regions of 
the cluster growth. From the applied point of view, stellated gold nanoparticles have nowa- 
days and projected in the near future, potential nanotechnological applications in disciplines 
such as the biomedical applications (cell tissue imaging, sensing and cancer therapy, drug 
delivery, etc.) 11-14 , and nanophotonics applications 41-44 . 
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FIG. 1: (a) Low resolution transmission electron microscopy image of gold nanoparticles with a 
high production rate of stellated particles. In (b) a magnified region where more detail of these 
stellated clusters is observed. 




FIG. 2: Transmission electron microscopy images of gold particles with assumed (a) cuboctahedral 
and (b) icosahedral seeds. Models that resemble these images are shown in (c) and (d). SEM image 
in (d) evidence both kind of stellated clusters coexisting. 



11 



C.L. Nehl, H. Liao, and J.H. Hafner, NanoLett. 6, 683 (2006). 



22 



0.5% 



0.4% 



0.3% 



HI 
-1 




0.2% 



0.1% 



B (10" 2 ) 



FIG. 3: Behavior of the fractional change in the cohesive energy accompanying a fractional change 
in the first nearest-neighbors distance e = (fij — ro)/ro as given by Eq. (4) for the potential models 
given in Table I. 
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FIG. 4: Some star-like gold clusters, similar to those shown in Fig. 2, where each one was obtained 
by using regular tetrahedral motifs attached to the eight (111) faces of a cuboctahedron (top row) 
or the twenty (111) faces of the icosahedron (bottom row), both with magic size. These pihata-like 
clusters correspond, from left to right, to k = 1, 2, 3, ..., 7, see Eqs. (6). 
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FIG. 5: Behavior of the ratio iV sur f/iV 2 ' 3 as function of N 1 ' 3 , where iV is the cluster size, for both 
the relaxed SC (A), and SI (Q) clusters, see Fig 4. The actual number of surface atoms of each 
stellated cluster, N mT f, is calculated through the Lee-Richard algorithm 35 using a probe radius 
r probe = 1-66 A. The dashed and solid lines are obtained by fitting the corresponding iV sur f data to 
Eq. (7). Note that the smaller cluster of the SI series is excluded from the fitting. 
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FIG. 6: Trends of the calculated excess energy per N 2 ' 3 , A(iV), as function of N ' 3 , for both SC 
(A), and SI (O) relaxed clusters for the potential models CR (a), CP (b), RGL (c), and BFFMM 
(d). The solid and dotted lines are the curves Asc(-^) 5 and Asi(iV) from Eq. (10), where the 
coefficients are given in Table III, respectively. For each potential model E re { = Asc has been used. 
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FIG. 7: Same as Fig. 4, but for the average equilibrium configurations of the stellated gold clusters 
obtained by MDS at room temperature. First row are those gold clusters based on initial structures 
with cuboctahedral core, and second row those with icosahedral core. The first three clusters of 
each series (k = 1,2,3) are structures of low-symmetry or amorphous-like clusters 34 ' 40 . The last 
cluster in each series corresponds to iVsi(9) = 6169 (bottom), and iVgc(10) = 5631 (top). 
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FIG. 8: Same as Fig. 5, but for the average equilibrium configurations shown in Fig. 7 and Fig. 11. 
In this case only the iV sur f data, for N Js 469 has been fitted to Eq. (7) excluding the coefficient 
7; the first three clusters of each series are excluded from the fitting, since they are amorphous or 
low-simmetry structures 34 ' 40 . iV sur f has also been computed through the Lee- Richard algorithm 35 
using a probe radius r pro b e = 1-66 A. 
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FIG. 9: Same as Fig. (6), but for the MDS results obtained using the CR potential model 19 ; the 
excess energy A(iV) has been calculated using the total energy. Note that we have excluded from 
the fitting the smallest clusters, since they are low-symmetry or amorphous structures with non 
star-like shape. 
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FIG. 10: The difference in excess energy or binding energy per N 2 ' 3 of the SC and SI clusters as 
function of cluster size. The curve is obtained from the fitted curves in Fig. 8. 
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FIG. 11: Snapshots of the equilibrated stellated cuboctahedral and icosahedral clusters. In (a) 
the clusters has iVsc(12) = 9437 atoms at T = 298 ± 0.46 K. The diameter of the cluster is 
around 10.8 nm; in (b) the cluster has iV S i(ll) = 10803 atoms. In this case T = 298 ± 0.44 K, 
and the diameter of the cluster is around 10.0 nm. To enhance surface details a molecular surface 
representation computed through the Lee-Richard algorithm 35 has been used. 
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